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We have obtained the solutions of two dimensional singular oscillator which is known as the quantum 
Calogero-Sutherland model both in cartesian and parabolic coordinates within the framework of 
quantum Hamilton Jacobi formalism. Solvability conditions and eigenfunctions are obtained by 
. . . using the singularity structures of quantum momentum functions under some conditions. New 
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^ i- In quantum mechanics, there are only hmited number of exactly solvable (ES) systems such that their eigenvalues 
and eigenfunctions can be obtained in an explicit and closed form [1,2,3,4,5]. Another definition is proposed that a 
quantum mechanical system is called as exactly solvable if the solutions of Schrodinger equation can be expressed in 

■ terms of hypergeometric functions mFn 0- There are some well-known potentials such as harmonic oscillator, Morse, 
trigonometric/modified Poschl- Teller, trigonometric/hyperbolic Manning- Rosen, and the Natanzon potentials 0,B0 

I— I ,, which are one-dimensional exactly solvable systems and their eigenfunctions can be stated in terms of hypergeometric 
i-C ' type functions i-Fi, 1F2 for bound states, 
^p^. On the other hand, quasi-exactly solvable (QES) problems have received great attention [10,11,12,13,14,15] in 
the last two decades. The separation of variables for the Hamilton-Jacobi and Schrodinger equations in multiple 
orthogonal coordinate systems is one of the chief problem in quantum mechanics, such as isotropic harmonic oscillator 
in three dimensions [16,17,18]. 

Q-i In the past few years, quantum Hamilton-Jacobi (QHJ) formalism [19j20,21,22] has been revisited. This formalism, 
formulated as a theory related to the classical transformation theory '2u\, was proposed by Leacock and Padgett in 
1983 and they also formulated QHJ method in one and three dimensions for a spherical potential which is separable 
[Tot . The application of QHJ to a class of solvable potentials has been re-explored in one dimensional space in 
' great detail by Bhalla, Kapoor and collaborators [2l|. It is shown that QHJ formalism is a useful method to obtain 
If^ ', eigenvalues and eigenfunctions of the Hamiltonian system, for both ES and QES problems [22]. Sextic oscillator and 
10 ' circular potential are discussed by Geojo and collaborators [12] as an application of the QES models. 

[ The motivation of this article is to introduce two dimensional QHJ approach as an adaptation of Leacock's work 
QQ i [1^ and to solve two dimensional singular oscillator which is introduced first time by Calogero-Sutherland [2^. In ref. 

■ [6], singular oscillator in parabolic coordinates is solved within polynomial expansion and solutions are obtained in a 
C , general form. In this study, as a first application of QHJ approach to singular oscillator, solutions will be obtained 

' both in cartesian and parabolic coordinates. Furthermore, we aim to generate new quasi exactly solvable potentials 
^ , by using Infeld-HuU factorization method [13] that has been developed recently by Dong [25j . 
V~j ■ Singular oscillator as the super-integrable potential is given as @ 

U' 21 

.5^; Viix,y)^^u;'{4x' + y') + hx+^2^ (1) 

where k2 > 0. It is stated in [6] that for fc2 > 1/2, the motion is restricted on one of the half planes (— 00 < a; < 
00, y > 0) or (—00 < X < cx), y < 0), whereas for < fc2 < 1/2 in whole plane (x, y). 

In the present work, in section II, two dimensional quantum Hamilton-Jacobi formalism is introduced and exact 
solutions are obtained for the singular oscillator in cartesian coordinates. Section HI involves QHJ formalism and 
QES solutions of singular oscillator in parabolic coordinates. In section IV Infeld-HuU factorization method has been 
used to generate new QES potentials. Conclusions are given in the last section. 

There are several formulations of quantum mechanics such as Heisenberg's matrix, Schrodinger wave function, 
Feynman's path integral, Wigner's Phase Space, density matrix, second quantization, variational, de Broglie-Bohm 
pilot wave and Hamilton Jacobi formulations. Quantum Hamilton-Jacobi formalism originally comes from the classical 
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Hamilton Jacobi theory which was first formulated in the form of present work by Leacock and Padgett in 1983 [l9|. 
Then this formulation was investigated in a different approach by Bhalla, Kapoor and their collaborators [21,22]. Let 
us write 

■p = VW (2) 
where W is the quantum Hamilton characteristic function related to the solution of the Schrodinger equation 

- 1^ VV(a:, y) + {V{x, y) - E) ^{x, y) = (3) 
by the definition i]j{x^y) = exp(^). Using this definition of wave function in Schrodinger equation gives 

iyWf - i^V .( V VF) = 2m {E - V{x, y)). (4) 
If Eq. ^ is used in Eq. (U the QH J equation is obtained as 

'p^ -ihV .~p :^2m{E~V{x,y)), (5) 
in other words 'p is the logarithmic derivative of ip{x,y) which can be written in the form of 

S/^p{x,y) 

p = 7 — r • 

v[x,y) 

Eq.® can be separated by using quantum Hamilton characteristic function W{x, y) = W{x, Ai) + W{y, A2). Here Ai 
and A2 correspond to two cartesian separation constants, E = Xi + A2. The quantum momentum functions (QMF) 
corresponding to the system are given as 

dW^ dWy 

from now onwards we set h — 2m — 1. Using Eq.([7]) in ([5]), QHJ equation can be separated as follows 

pI - idxPx - Ai + 2cj^.T^ + kix = Q (8) 



pI - ^dyPy - A2 + y y2 ^ ^ . (9) 



To solve Eq.® within the QHJ formalism, let us change the variable as z = x + Hence, Eq.® turns into the 
form of 



pl{z)-^^-{e^-2u;^^J)^0 (10) 

where ei = Ai + . In Eq. pTI)) . QMF pz{z,Si) = Pz{z) has n poles which are the zeros of the wave function. As 
it is seen from the Eg. pU)) . the residue at each of these poles is — i. This QMF has not any other poles except at 
infinity. In the limit of large z, QMF approaches to Pz{z, ei) « ±i\/2ujz. Then we can write 

Pz{z,ei) K ±iV2ujz + Q{z) (11) 

where Q{z) coming from the Liouville theorem is to be determined later. The correct sign oi pz{z, ei) can be obtained 
by using the condition of square integrability of the wave function. The wave function is given as ipiz) = f P^i^^'^i)d-z 
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for a bounded wave function at large z, we choose the positive sign of QMF in Eq. ([TT|l . Then one can write the QMF 
function as 



E 

fe=i 



+ i\/2ujz + (t>{z), 



(12) 



where (j){z) is analytic and bounded at infinity. From Liouville's theorem [19,20,21,22], (l){z) has to be a constant 
which 
yields 



which is zero for z -> C50. In Eq.([I2D ^ ^ = where P{z) = JJ {z — Zk). Substitution of Eq. fT^ in Eq.lT 



k=l 



k=l 



P"{z) - 2V2ujzP'{z) + (£i - V2lo)P{z) = 



(13) 



the solution of this equation gives the Hermite polynomials Pn{z) = H„{\/V2ujz) and Ai can be obtained as Ai = 
{2n + l)V2uj- The wave function has the following form 



ipiz) = exp(i J pz dz) 

V2 I 

= Bn exp[ — ^cjz^] iJ„[Y V2a;z]. 

where -B„ is a normalization constant. Now let us consider the solutions of Eq.(9); there is a singularity at ?/ = in 
Eq.®, then, one can expand = ^ + oq + aiy + 022/^ + ... and use in Eq.Q, therefore following quantities can be 

VI. 1 1,2 , 1 



(14) 



found: oi = , = and 61 = — ± 2^ y ^2 + i- If we expand Py as 

bi iuo .P „ 



where P is an nth degree polynomial, prime denotes the derivative in terms of y and C appears from Liouville theorem 
[19,20,21,22] which equals to zero for large y. By using Eq. fMl) in Eq.®, the following form can be obtained: 



So we get the solutions of this equation as 

A2 = t^\/2 

and 



= 0. 



(2n+l)T^\/fc: 



1 



(15) 



(16) 



Pn{v) 



(zb, - i) 



T ^' 



(ibl-l) / ^ ,2 



2 



V2 



y 



where L'^(y) stands for the associate Laguerre functions and (a)„ :— a{a + I). ..{a + n — 1) for n > 0, 
(a)o = 1 [24|. Using Pniy), the wave function is obtained as 



here A'^ is a normalization constant. Using 



in Eq. pTp . then wave function becomes 



iV2 /T 



My)^N 



2V2 



y')L 



V2^ 



(17) 



(18) 



(19) 
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When bi has a positive sign, then k2 < ^. As a special case; for the value of ^2 = 5, positive sign of bi should be 
taken and in this case Laguerre functions can be written in terms of Hermite polynomials by using the relations [2^ 

H2n{y) = (-l)"22"n!i,-i/2(y2) ^ H2n+i{y) - (- l)"22"+in!yi,- (y2 ) . (20) 
From the values of Ai and A2 , the total energy of the system can be written as 



£;„ = V2uj{2n + 2) - ^ T fci + 4 

where Ai and A2 are labeled with rii and 71,2, n = ni + 712. The results agree with the difference of arises from 
choosing the separation constants. 

Parabolic coordinates ^ and rj are connected with the cartesian x and y hy x — ^(^^ — ?7^) and y — ^tj respectively, 
for ^ g R, r/ < 0. Then, the potential relation in Eq.([T]) turns into the following form 



Schrodinger equation in parabolic coordinates is given as 

- 1^ vV(^, V) + [vit v) - E) ^{t V) = 

2 m 

where = g2^^2 + By using the definition tp{£,,r]) — exp(^) in Schrodinger equation gives 

(VWf - ihV.(yW) ^2m{E- ^(^,77)). 

The logarithmic derivative of "0(^1^) gives the quantum momentum function ~p* = — From now on we set 
h = 2m = 1. Then, quantum characteristic function T4^($, 77) and quantum momentum functions are defined below 

Wit 77) = Wit E,) + Win, E2), Pi = p, = ^ 

Hence, QHJ equation in parabolic coordinates can be separated as 

pI - iPi - E^e + + ^fci^^ + ^-T = Q (22) 

pI - - E2if + yT]^ - \kiri^ + ^+T = (23) 
where A = ^^-^'^ g^j^^j T is a separation constant. Let us solve Eq. (P^ and firstly expand in series as 

p(e,^i) - j+ao + aie + a2e' + a3f (24) 

and substitute in Eq. ([22|) . we obtain 03 = ±i^/B, 02 = 0, ai = ±-^^, qq = zt-v/T and b\ = 1 + e^)- Here, 
B — £ = ±1 and /i = \f\\~AA. Now we expand in Laurent series as 

p(e,i^l) = ^±VT+-%^ + ^^/Bf -*§i + Ci (25) 

where Pnit) is an 7T,th degree polynomial which is unknown and C\ is a constant due to Liouville theorem. If we use 
Eq. ([25| in Eq. ([22|) . for large values of ^, P = ^" and this leads to Ci = 0, T = 0. This case also leads to 



E, ^ ^u^in + 2) ± i=^i + fc| - ^ (26) 
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in ^ — > cx) limit. The main equation is obtained by using Ea. (|25p in Eq. ([22l 



8^ 



2VB 



16B 



2ibiVB-Ei]^ 



ki ibiki 



PAf(e') = 0. (27) 



where Pm(C^) denotes the polynomials of degree M. With this form, this equation can not be solved conventionally 
given in [22]. The difference in our solution arises from the parameter transformations. We use 



El = 'sVb ^ + 2ibiVB + aVbM 

16B 



in Eq. ((27|) for notational reasons. Thus, this equation becomes 



+ 2VBe -'^]d^~ AVBAie - e + -^(1 + 3z6i) 

2Vb ^ J ^ aVb 



(28) 



(29) 



where e = ^J^- The simplest solutions of Eg. ([29]) will be introduced. 

Case 1: Now let us find one explicit solution which occurs when A/ = 0. In this case, the polynomial 

^o(^') = 1 

is a constant. Then, 

+ fci / 5 



2uj \2V2 4 



(30) 



(31) 



is obtained for the positive sign of the residue hi. This also leads to ^2 = Corresponding wave function is given 
as 



+ 



(32) 



Here, the boundary conditions require ^/l + 4fc| < i, but using fc2 = -3^7 it is seen that the solutions arc not 
physical for the wave function. And the solutions for the negative sign of the residue bi are 



fci / 5 



2lu \2V2 4 



that leads to fc2 = ^ and 



(33) 



(34) 



Case 2: AI ~ 1 stands for the solutions in two dimensional space which requires 



where s is a constant. Substituting this into Eq. ()29|l . one obtains 

3 



/ k^ 
eia = ^^(1 + T Y 8/B + l&ibi^/B + ^ 



(35) 



(36) 



Again we can consider the solutions according to signs of the residue. For the positive sign of residue, solutions are 
obtained as 



1 



2\/2w V2 2V2 



1 + 4fc2 T 7 



(37) 



(38) 
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where 7 — J %y/2L0 — Auj^Jl + 4/c| + ^jj- and — y= + 4fc| < i. Solutions for the negative sign of the residue are 



2^ 2\/2 V ' ^"-2 ^ 2 

introduced as 

3 /3 1 



^- = ^^+57fVi+4fci J^l^fc^ ± ^^^^ ^2 _ 2x/2tj (^2 + -^^^1 + 4^1^^ e"(^*'+^^') (40) 
where 7 = ^^8^2^ + 4^^1 + 4/1;^+^ and i + ^VlTIfcf > 0. 



One can transform Ea. (|29p into a form 



by using 



/"(f) + f f 3x/B - ^ - 2^b,V^ + 4x/BA/) f - Be ~ + e ] /(C^) = (41) 



P(f ) = e^—^ )/(f ) (42) 



where e = — e + "^^7^ + O^^c can see that Ea. (|¥T|l is in the form of Schrodinger equation with a potential 

V{0 = Bf + - (sVs - ^ - 2z6i + aVbM^ e (43) 

For large values of M, solutions of Eq.(|4T]) can be obtained but an ansatz for the / which is more compact can be 
written as 

i=l 



where ^/rji are the zeros of the /(f) and <; — ^^^3/2 ~ i^i + f ■ Then, e is obtained by using Ea. (|l^ in Ea. (|¥T|) as 

fci _^ ^ 
8V2w 2 



32B3 

A/ 



e' = (2Af + 1)a/2^ + ± + 4A:| - 4^2^; ^ r?,. (45) 



Following the same procedure, the solutions of Ea. ((23)l can be obtained as in the previous way. 

Now we introduce Infeld-HuU factorization method firstly [25,26]. For a given QES V{S) potential, Schrodinger 
equation is [15] 

^'■{0 = {V{0-E,)^,. (46) 
Using following logarithmic derivative in Schrodinger equation 

gives Riccati type equation 

V{0-Ej^y^{0+y^,{0- (48) 

For M + 1 solutions Schrodinger equation is 

^2fc+, (0 = {V{i) - E2k+M2k+, (49) 

/c = 0, 1, ...M [15] and from Eq.®, Eq.(|171) and Eq.® we get 
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If the operator — ^ acts on 
then, the Schrodinger equation with a new potential and wave function is obtained [15]: 



(yj + l^^iyj ~ 77^)V'2fe+j = {E2k+j - Ej)'iP2k+j (51) 



(-^ + ym^2k+,{o = E2k+,^J2k+m (52) 



where 



V{0 = V{^)-2^lni,,{^) (53) 



7 (T? ^ jioo^2k+j^,d^ 

^ = {E2k+j - Ej ) — . (54) 

From Eq.(53) it seems that new potential is generated by Darboux type transformation from V(£,) [27,28,29]. In 
the language of supersymmetric quantum mechanics V^(^) and V{^) are known as supersymmctric partner potentials. 
Now new potentials can be obtained for M ~ and M — 1 cases using Eg. ((53|) . For the case of M ~ 0, new potential 
is 

V = + l-E + 3V2c. + ^)e + ^ ^ ^ + ^ ^ (55) 

where ± arises from positive and negative signs of the residue values. For M = 1, the potential has a more complex 
form which is 



where 



^ 2(&V:B + ^)^ 4a^(2+p)V B 
2 ^ «2+6)2 (0^2+5)2 



^6(0 = — + /„.2^M2 r.^2^M2 ' (57) 



WS(2^+7)-i-^ 4abp + 4i^+bVB)iai2+p)^^) 

^^(^)^T + <^T6 WTW ' ^ ^ 

(-^^ + ||-2V^bp(2+|)) 2a^(2+p--liL)^-2(5^pV^+^ + 5^pV^) 

'^^^^) = -^ <^T6 + WTW ' ^''^ 

= ^ <^+K^' ^ ^ 

2a(2+rt(l+p)-^ 
= , (61) 

.02(0^^^^.^^, (62) 
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where a = ki±ujV2, b = -2\/2w(2± -^y/TTIk^), P=\± 

As a first application of two dimensional QHJ approach, we investigate the solutions of 2D singular oscillator 
potential which admits separation of variables in two different coordinate systems. After separation of variable, this 
potential can be handled as two different one-dimensional systems. Then, by following one dimensional QHJ approach 
of Bhalla and collaborators work, we get the exact solutions in terms of Hermite and Laguerre functions for cartesian 
coordinates which are consistent with the solutions in [s^l ■ In section III, it is shown that the present method [l^l is 
problematic for obtaining polynomial solutions in parabolic coordinates. Hence, we use transformations of parameters 
in order to obtain solutions for M = 0, 1 cases, thus only the limited number of solutions are obtained in parabolic 
coordinates. The transformation of Ea. (|29p leads to a Hamiltonian with a new potential and solutions of this form are 
obtained in more compact form. Despite the more general results obtained in Ref.[6] by using polynomial solutions, 
it is stated that a limited number of wave function solutions can be obtained correctly by using QHJ approach for 
singular oscillator in parabolic coordinates. Thereof, it can be an open problem to investigate a generalization of the 
approach for the QES problems. 

Finally, new quasi-exactly solvable potentials are constructed for M = and M — 1 cases by using Infeld-HuU 
factorization method. By means of this method the solutions of more complicated potential can be found just by using 
the solution of singular oscillator potential from Eq.(54), so there is no need to solve the corresponding Schrodinger 
or quantum Hamilton- Jacobi equation. 

This work was supported in part by the Scientific and Technical Research Council of Turkey (TUBITAK). 
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